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Abstract 

This is a set of four lectures devoted to simple ideas about turbulent transport, a ubiquitous non- 
equilibrium phenomenon. In the course similar to that given by the author in 2006 in Warwick [45], 
we discuss lessons which have been learned from naive models of turbulent mixing that employ simple 
random velocity ensembles and study related stochastic processes. In the first lecture, after a brief re- 
minder of the turbulence phenomenology, we describe the intimate relation between the passive advection 
of particles and fields by hydrodynamical flows. The second lecture is devoted to some useful tools of the 
multiplicative ergodic theory for random dynamical systems. In the third lecture, we apply these tools to 
the example of particle flows in the Kraichnan ensemble of smooth velocities that mimics turbulence at 
intermediate Reynolds numbers. In the fourth lecture, we extends the discussion of particle flows to the 
case of non-smooth Kraichnan velocities that model fully developed turbulence. We stress the unconven- 
tional aspects of particle flows that appear in this regime and lead to phase transitions in the presence of 
compressibility. The intermittency of scalar fields advected by fully turbulent velocities and the scenario 
linking it to hidden statistical conservation laws of multi-particle flows are briefly explained. 



1 Turbulence and turbulent transport 

Hydrodynamical turbulence, with its involved flow patterns of interwoven webs of eddies changing errati- 
cally in time [68] , is a fascinating phenomenon occurring in nature from microscopic scales to astronomical 
ones and exposed to our scrutiny in everyday experience. Despite the ubiquity of turbulence, an under- 
standing of this far-from-equilibrium phenomenon from first principles is a long-standing open problem, 
one of the few remaining theoretical challenges of classical physics. Although it seems still a too far- 
fetched goal, there has been a constant progress over years in the theoretical and practical knowledge of 
hydrodynamical turbulence and, even more so, of other non-equilibrium phenomena of a similar nature. 
This progress has been due to developments of experimental techniques and computer power, but also 
to new theoretical ideas inspired by studies of simple models of turbulence-related systems. The present 
lectures will discuss one circle of such ideas more relevant to the problem of transport properties in tur- 
bulent fiows than to turbulence itself. They concern properties of stochastic processes underlying the 
turbulent mixing, see also ,75. 39j. In refs. [601 1781 [751 1791 176| one may find discussion of other theoretical, 
experimental, numerical and practical aspects of turbulent transport. For the use of stochastic processes 
in modeling of non-linear fluid dynamics itself, see e.g. [32]. 

It is believed that realistic turbulent flows are described by hydrodynamical equations, like the Navier- 
Stokes ones, that govern the evolution of the velocity fleld and, eventually, of other relevant fields like 
fiuid density, temperature, etc. As any structure whose dynamics is governed by evolution equations, 
the turbulent flow may be viewed from a theoretical point as a dynamical system. Many attempts were 
made then to apply to turbulence ideas from the dynamical systems theory, see e.g. [22| . These ideas, 
developed for low-dimensional dynamical systems, although often successful in explaining the onset of 
turbulence, seem to miss important features when confronted with the fully developed turbulence. Some 
of them, however, prove very useful in the description of transport in turbulent fiows, as we shall see in 
the sequel. 
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1.1 Navier-Stokes equations 



These are equations dating back to the work of Claude-Louis Navier from 1823 and of George Gabriel 
Stokes from 1883, completing Leonhard Euler's equations of hydrodynamics from 251 years ago. They 
have the form of the Newton equation for the mass times acceleration of the the fluid element: 

p{dtv + {v ■ 'V)v) — pi^'V'^v — —'Vp+f . 

/ T / \ 

fluid kincmatical pressure foree 

density viscosity density 

and have to be supplemented with the continuity equation 

ftp + V ■ {pv) = 

and an equation of state F{p,p) = 0. For the incompressible fluid, the latter states that the density 
p = po is constant and the Navier-Stokes equations may be rewritten in the form 

dtv + {v-'V)v - uV^v = ±{-Vp+f) (1.1) 

and are supplemented with the incompressibility condition V • i) = 0. In the latter case, they may be 
viewed as an evolution equation on the (infinite dimensional) space ^ of the divergence-free vector flelds 
describing an infinite-dimensional (non-autonomous) dynamical system 

^ = ^{t,X) forXG^o. (1.2) 

In the theoretical modeling of the way in which the fluid motion is induced, one often assumes that 
the external force / is random. The right hand side ^'{t,X) will then be also random, describing an 
infinite-dimensional random dynamical system. 

Eq. (|l.ip is nonlinear. The strength of the non-linear advection term (d • V)t; relative to the linear one 
u'V^v describing the viscous dissipation depends on the length scale I and is measured by the running 
Reynolds number deflned as 

V 

where A;u is the size of typical velocity difference across the distance scale with r; = -^-^ giving the 
typical turnover time of eddies of size I. In particular, if L is the integral scale corresponding to the 
size of the flow recipient, Re{L) = Re is the integral scale Reynolds number and tl is the integral time. 
On the other end, Re{rf) = 1 for r] equal to the viscous (or Kolmogorov) scale on which the non- linear 
and the dissipative term in the equation have comparable strength. The Kolmogorov time gives the 
turnover time of the smallest eddies. Phenomenological observations of hydrodynamical flows point to 
the following classification according to the size of the Reynolds number Re: 

• ^6^1- laminar flows, 

• _Re ~ 10 to 10^: onset of turbulence, 

• > 10^: developed turbulence. 

In he laminar regime, some explicit solutions are known |1U) . At the onset of turbulence, the flow is 
driven by few unstable modes. The standard dynamical system theory studying temporal evolution of 
few degrees of freedom governed by ordinary differential equations or map iteration proved useful here 
e.g. in describing the scenarios of the appearance of chaotic motions, see e.g. [3S]. Finally, in the fully 
developed turbulence, there are many unstable degrees of freedom. Kolmogorov's scaling theory [54] for 
this regime predicts that 

Aiu oc l^^^ for -q^l^L 

which is not very far from the observed behavior. The number of unstable modes may be estimated 
as given by {L/rff oc Re^^^. New phenomena arise here that have to be addressed, like cascades with 
(approximately) constant energy flux, intermittency, etc. 

1.2 Turbulent transport of particles and fields 

One may view a turbulent flow as a dynamical system in another way, related to the transport of particles 
and fields by the fiuid 
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1.2.1 Transport of particles 

The particles may be idealized fluid elements, called Lagrangian particles, or particles without inertia 
suspended in the fluid that also undergo molecular diffusion. We talk in the latter case of the 

• tracer particles 

whose motion is governed by the evolution equation 

d ?? 

^=v{t,R) + V2;^rj{t), (1.3) 

where k stands for the (molecular) diffusivity and r]{t) is the standard vector- valued white 
noise. The Lagrangian particles correspond to the case k — 0. 

Small objects like water droplets in the air experience a friction force proportional to their relative velocity 
with respect to the fluid [61) . These are 

• particles with inertia 

whose position R and velocity U satisfy in the simple case the equations of motion 

^ = ^ = L{^U + v{t,R) + V2;^ri(t)), (1.4) 

where t, the Stokes time, accounts for time delay with respect to tracer particles whose 
evolution is recovered in the r ^ limit. 

Finally, one may consider small objects with spatial structure suspended in the fluid, like 

• polymer molecules 

whose motion is sometimes modeled [20] by the differential equations 

^=v{t,R) + V2i^r]{t), ^ = (S ■ V)v{t, R)-aB + V2^i{t) , (1.5) 

where R describes the position of one end of the polymer chain, B is the end-to-end sep- 
aration vector, and ri{t) and ^{t) are independent white noises. The term ~aB describes 
the elastic counteraction to the stretching of the polymer. 

All three cases provide examples of flnite-dimensional dynamical systems if w,t7,^ are given. The dynam- 
ical systems are random if d, t;, ^ are random with given statistics. They may be studied with tools of the 
dynamical systems theory to which they have also provided new inputs. 

1.2.2 Transport of fields 

Turbulent flows may also transport physical quantities described by flelds changing continuously from 
point to point, like temperature, pollutant or dye concentration, magnetic fleld, etc. The evolution of 
such transported flelds is described for small field intensity by linear advection-diffusion equations with 
variable coefficients related to the fluid velocity. Under the same assumption of small fleld intensity, one 
may ignore the influence of the transported fleld on the fluid dynamics, ending up with a passive transport 
approximation. The passive advection of a 

• scalar field 9{t,r) 

(e.g. temperature) is governed by the equation 

dte + {v-'V)e-K'v^e = g (i.e) 

where k stands for the diffusivity of the scalar and g is a scalar source. 

Similarly, for a 

• density field n{t, r) 

(e.g. of a pollutant), one has the partial differential equation 

dtn + V ■ (nv) - kV^u = h (1.7) 
that differs from that for the scalar field only for compressible velocities with V ■ « 7^ 0. 
For the passive transport of a 
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• phase-space density n{t,r,u) 

(e.g. of an aerosol suspension of inertial particles) one can write the equality 

dtn + (u • Vr)n + V„ ■ (n " ~ ^ ) - kV^ti = h . (1.8) 

r 

Finally, the passive transport of the 

• magnetic field B(t,i) 

satisfying V ■ B = is governed by the evolution equation 

dtB + {v ■ V)B + (V ■ v)B - {B ■ V)t; - kV^B = G , (1.9) 

where k is the magnetic diffusivity and G is a source term. 

The meaning of "small intensity" in the assumptions leading to the above equation depends on the 
physical situation (e.g. very small admixtures of polymers or chemical reactants may change the flow 
considerably) . 

1.2.3 Relations between transport of particles and fields 

The dynamics of localized objects (particles, polymers) and passively transported fields listed above are 
closely related. 

Let Ji(t; to, ^o|»7) denote the time t position of the tracer particle that at time to passes through ro. 
It depends on the realization of the noise r], see Eq. (|1.3p . Then the formula 



e{t, r) = e{to, R{to; t, r jj?)) = J e{to,ro) 5{ro - -R(to; t, r\r))) dro , (1.10) 

where the overline stands for the average over the white noise 77 realizations, expresses the solution of 
the advection-diffusion equation (|1.6[) with vanishing source g in terms of the initial values 6{to,r) of 
the field. The meaning of this solution is that scalar field is conserved along the noisy trajectories. 

Let us consider the matrix W{t\ to, ro\r])) that propagates the separations between two infinitesimally 
close tracer particle trajectories: 

W j{t;to,ro\'n) = —-J ■ 

The determinant of matrix W describes the flow-induced change of the volume element along the tracer 
particle trajectory and it is equal to 1 for the incompressible flows. The formula 



n{t,r) = AetW(to;t,r\'n) n{to,R{to;t,r\-n)) = j n{to,ro) S{r - R{t;to,ro\Ti)) dro (1-11) 

gives the solution of the initial-value problem for the density evolution equation (|1.7p with vanishing 
source h. Note that 

W{to;t,r\ri) = W{t;to, R{to;t,r)\riy^ 

so that the density changes inversely proportionally to the change of volume element along the Lagrangian 
trajectories. The density increases where the flow contracts the volume. 

Similarly, denoting by {R{t; to, ro, uo\ri), U(t; to, ro uo\ri)) the phase-space trajectory of the inertial 
particle passing at time to through (ro,Uo) and introducing the matrix W{t;to,ro,uo\ri) propagating 
the phase-space separations between two infinitesimally close inertial particle trajectories, one obtains the 
formula 

n{t,r,u) = det W {to; t, r, u\r]) n{to,R{to;t,r, u\r]),U {to; t,r,u\r])) (1-12) 

for the solution of the evolution equation (|1.8|) with h = 0. 

Finally, the solution of the evolution equation (|1.9p for the magnetic field with G = is given by the 
formula: 

B{t,r) = {detW{to;t,r\Tj)) W{to;t,r\ri)-^ B{to,R{to;t,r\Ti)) . (1.13) 

The action of the matrix W{to;t,r\ri)~^ = W{t;to, R{to;t,r\T])\T]) stretches, contracts and/or rotates 
the magnetic field vector along trajectories. Note in passing that the solution of the second of Eqs. (II. 5|) 
for the polymer end-to-end separation in the linear approximation is given by 

B{t) ^ e-"^'-'°^W{t;to,ro\'n)B{to) (1.14) 
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if the component a of the polymer diffusivity vanishes. 

The source terms providing inhomogeneous terms in the linear advection-diffusion equations for the 
fields may be taken into account by the standard "variation of constants" method. For example the 
solution (|1.10|) of the scalar evolution equation (II. 6|) picks up for non-zero source g an additional term 



on the right hand side. 

Conclusion. Passive transport of particles and fields by hydrodynamical flows, including turbulent ones, 
IS described by simple equations closely related m both cases. The transport of particles will be studied 
below by the (random) dynamical system methods. We shall subsequently examine which attributes of the 
particle dynamics bear on which properties of the field advection, with the goal to capture and explain 
some essential features of the turbulent transport. 




(1.15) 
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2 Multiplicative ergodic theory 



We shall be interested in a statistical description of particle dynamics in turbulent flows. For concreteness, 
we shall concentrate on the dynamics of Lagrangian particles carried by the a flow in a bounded region 
V. Their trajectories are solutions of the ODE 

^^v(t,R). (2.1) 

with random, not necessarily incompressible, velocity fields v{t,r), see Eq. (|1.3p . We shall assume that 
the statistics of the velocities is stationary. The deterministic time-independent velocity field will be 
also covered as a special case. Since Eq. (|2.1|l has a form of a general dynamical system, most of the 
considerations that follow will apply without much change also to other particle-transport equations. For 
reasons which were already mentioned in the introduction and will become clear in later, the considerations 
of the present and the next lecture, borrowing on the theory of differentiable dynamical systems, are 
relevant to the advection by turbulent flows at intermediate Reynolds numbers. 



2.1 Natural invariant measures 

In deterministic time-independent flows, one calls a measure n(dr) on the flow region V invariant if 

j f{R{t;r)n{Ar) ^ j f{r)n{dr) (2.2) 

V V 

for bounded (measurable) functions f on V and all times t, with the shortened notation R{t\ r) = 
R{t; 0, r). This may be generalized to the case of stationary random flows by considering a collection of 
measures n{dr\v) on the flow region, parametrized by the velocity realizations, such that 

J f{R{t;r\v)\vt) n(dr|t;)^ " (/ '^^'^^'"^ "-(dr|t;)^ = J f{r\v) 7V(dr|dt;) (2.3) 

V V VxV 

where ( — ) denotes the average over the velocity ensemble and Vt{s, r) — v{s ~\- 1, r). Such a collection 
defines a measure N{dr\dv) on the product V x V of the flow region V and the space V of velocity 
realizations. This measure is invariant under the 1-parameter group of transformations 

VxV 3 {r\v) I — > {R{t-r\v),v) . 

It is called an invariant measure of the random dynamical system 

Let us seed the particles into the fluid at some past time to with smooth normalized density nto {to, ro), 
e.g. the constant one. This density will evolve in time according to the advection equation dtn-\-'V -{nv) = 
0, see Eq. (|1.7p . In an incompressible flow, the constant density nto{t,r) = IV'I"^ will stay constant and 
normalized also at all later times t. If the flow is compressible, however, then the particles will develop 
preferential concentrations and the densities nt„{t,r\v) will become rougher and rougher with time in 
typical velocity realizations as illustrated by the snapshots taken form [21) : 



Similar phenomenon occurs for the phase-space density of inertial particles, see e.g. [121 113j . The proba- 
bility measures n{dr\v) on V are called natural measures if 

f{r)n{t,dr\v)^ lim ttt/ //(»") "to (0,dr|w) 

fo^-oo to J. J 
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for all continuous functions / and almost all velocity realizations v. In particular, / f{r) n{dr\v) = 
lim J /(r) ntg (0, r\v) dr if the latter limits exist. In that case, the natural measures are just the time 

to — » — OO 

zero distributions of particles seeded with a smooth density at time — oo. For compressible flows, the 
natural measures n{dr\v) are concentrated on a random time- and velocity-field-dependent attractor. 
These measures have the property l|2.3p and give rise to the natural invariant measure A'^(dr|d?;). The 
latter permits to do statistics of the values of functions f{r\v) of points in the flow region and velocity 
realizations. Such statistics is often called "Lagrangian" since it may be calculated from time-averages 
along the Lagrangian trajectories. 



2.2 Tangent flow 

Recall that the matrix with the elements 

_ dRHt;to,r\v) 



propagates infinitesimal separations 5R between Lagrangian trajectories in a given velocity realization: 

SR{t-to,r\v) = W{t-to,r\v)Sr . 

We shall abbreviate W{t;0,r\v) = W{t;r\v) dropping also the other dependences in W from the 
notation whenever they are not necessary. Diagonalizing the matrix W^W and writing its positive 
eigenvalues in the exponential form as e'^''^ > e^''^ > ■ • ■ > e'^'''*, we obtain the so called "stretching 
exponents" pi{t;r\v) with the ordering 

pi > p2 > ■■■ > Pd- 

Due to the cumulative effect of short-time stretchings, contractions and rotations of the separation vectors 
of infinitesimally closed Lagrangian trajectories, the stretching exponents typically grow in absolute value 
in time. We shall be interested in their long-time asymptotics. Let us call the ratios ai — the 
"stretching rates". They are sometimes also called "finite-time Lyapunov exponents". For fixed time t, 
the functions ai{t;r\v) may be treated as random variables on the product space x V equipped with 
the natural invariant measure A'^(dr|dt;). Their joint probability density function (PDF) is then given 
by the formula; 

Pt{a)= J 5{a~a{t;r,v))N{dr\dv). 

vx-r 

The PDF's Pt{S) and P^tiS) are very simply related: 

Pt(a) = P-t(-5), (2.4) 

where a = ((t^, . . . , tJi) if <? = (cri, . . . , a^). This follows from the relation R[—t; R{t; r\v)\vt) — r 
implying by the chain rule that W{t;r\v) = W {~t; R{t,r\v)\vt)~^ and that 

a{t;r\v) = ~a{-t- R(t,r\v)\vt) (2.5) 

and from the invariance (|2.3p of the measure N(dr\dv). A very general result, the backbone of the 
multiplicative ergodic theory, assures that the stretching rates approach constant limiting values for long 
times. Denote ln+ x = max(0, In a;). 

• Multiplicative Ergodic Theorem (Oseledec 1968 [66], Ruelle 1979 [70]). 

If the measure invariant measure A'(dr|d?;) is ergodic and /ln+ || W^(t; r |?;) || A*'(dr, dw) < oo for 
some t > then 

lim Pt{a) = 5{a~X) (2.6) 

t — >oo 

for a vector A — (Ai, . . . , Ad) of "Lyapunov exponents" such that Ai > • ■ ■ > Aj. 

This may be viewed as a particular case of a multiplicative version of the law of large numbers [5]- For 
typical r, v and Sr, 

Ai = lim - ln|lVK(t;r|-i;)5r|| . 

t — 'OO t 

The strict positivity of Ai signals the sensitive dependence on the initial conditions and is taken as a 
definition of "chaos". Similarly, for 1 < n < d and typical Sri, . . . , 5r„, 

Ai + ---+A„= lim ^lndet(5r,-W'^W{t;r\v)5rj)i<,,,<n). (2.7) 

t^oo At \ / 



It is important to know how the limit (|2.6|l is approached. If Ai > ■ • ■ > A^ then one may expect that 
the fluctuations of the stretching rates about their limiting values behave according to 
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Multiplicative central limit (not covered by a general theorem) asserting that 

lim ri/^p,(A + r^/V) 



t^oo ^ ^ det(27rC)i/2 ■ 

Finally, as argued in [7], see also [31] , again under the assumption that Ai > • ■ ■ > Ad, one may expect 
occurrence of the regime of 

• Multiplicative large deviations (again not covered by a general theorem), where 

- lim i In Pi (a) = H{a) 

for a "rate function" H > such that H{X) = 0. 

For incompressible velocities where ai = 0, the rate function H{a) is infinite unless a satisfies the 
latter condition. Loosely speaking, the existence of the regime of large deviations means that 



Pt{a) « e 



-tH{S) 



for large t. As in the additive case, the existence of such regime with the rate function H{(t) twice 
different iable around the minimum implies the central limit behavior with 

q2 rr 

doiduj 

There are partial results about existence of the multiplicative large deviation regime for deterministic 
dynamical systems [SD] and for random dynamical systems with decorrelated velocities [111 131] . The 
latter case will be discussed in the next lecture. The large deviation regime for the stretching rates in 
realistic flows becomes accessible numerically and even in experiments [231 [9]. 

Let us denote by v' the time-reversed velocity field: v'{t,r) = —v{t,r). The large deviations rate 
function H(B) and H'{a), the latter one pertaining to the fiow in the time-reversed velocities, should 
satisfy 

• Multiplicative fluctuation relation (of Gallavotti-Cohen type) 

d 

H'{^a)=H{S)-Y,ai (2.8) 

i=l 

In particular, for time-reversible velocities when the fields v and v' have the same distribution, the two 
rate function coincide: H' = H. The original fiuctuation relations studied the large deviation regime of 
s = — X^f^i pi equal to the (phase-)space contraction exponent. The latter has been interpreted [711 172] 
as the entropy production of the dynamical system, with a = s/t standing for the entropy production 
rate. For deterministic uniformly hyperbolic (discrete-time) systems that are time-reversible, Gallavotti 
and Cohen [33] established that the rate function /i(ct) of large deviations of s satisfies the relation 

h[-a) = h{a) + a . 

This was generalized recently in [251 to random uniformly hyperbolic systems. Since 

h{a) — min H{a) , 



the relation (]2.8[) implies Eq. (]2.2p . 

Let us prove here following [6], see also [39], another, simpler, relation for a modified joint PDF of the 
time t stretching rates defined by 

Pt{ff) = (^ls(a-a{t;r\v))^'j (2.9) 

Note that Pt employs to average over r rather than the natural measures n{dr\v) used in Pt. 
Numerically, it is in general simpler to attain the PDF Pt than Pt as the latter requires performing 
longer Lagrangian averages. 
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• Transient multiplicative fluctuation relation (of the Evans-Searls type) 

P;(-a) = A(a) (2.10) 

where P' (S) pertains to the time-reversed flow. 

The proof of the relation (|2.10p employs a simple change-of-variables argument similar to the ones used 
to obtain transient fluctuation relations [361 1371 [53] . First, using the definition (|2.9p and the stationarity 
of the velocity ensemble, we obtain 

P_,(-a) = l^j 5{-a-a{-t-R\v))^^ = (^J 5(a + 5(-t; ^ ) • (2.11) 

Upon the substitution R = R{t; r\v) with the Jacobian 



this gives 



^ = detW{t;r\v) = e''! "-'(^^-W (2.12) 



y5(a-a(t;rW)el'ISS'^^'^'(*^'-|'')^^ = P,(a) e'*! ^^-^ 



where we have use the relation (|2.5|l to obtain the second equality. Now, the straightforward relation 
R{—t;r\v) = R{t;r\v') implies that 

P„,(a) = P;(a). (2.13) 

and the identity (|2.10|l follows. 

For long positive times t ^ 0, the Lagrangian particles approach the (time-varying) attractor and 
their statistics should not depend much on whether their initial points were distributed with the uniform 
measure or with the attractor measure. One may then expect that 

Pt{a) ^ Pt{a) (2.14) 

at long positive times, and similarly for P/ and P/. If this holds and the PDF's Pt and P/ exhibit large 
deviations regime then (I2.10|) implies (|2.8|) . Both assumptions are not granted, however, and require a 
careful analysis in specific situations. 

2.3 Applications of multiplicative large deviations 

Although concerning relatively rare events, large deviations play an important role in turbulent transport 
phenomena. Indeed, rare fluctuations of concentrations of toxic emissions or chemical agents may deter- 
mine the nocivity and/or the reaction rates. Suppose that we put at time zero a blob of density no(j'o) 
into a compressible fluid. Then, for vanishing diffusivitj0, the blob density at later times is 

n{t,r) = detW{0;t,r) no{R{0;tr) , (2.15) 

see Eq. (|l.ll|l . Along the trajectory R{t;ro), 

n{t) = n{t,R{t;ro)) = det VF(f; ro)" Vo(ro) = e" ''•'*'''o'no(ro) 

and 

{n{ty) = noiro)" j e""* ^ Pt(a) da « no(ro)' j e""* ^ da oc e"^"* (2.16) 
for large t, where 

-7,= min qy^ai + H{a). (2.17) 

As we see, the asymptotic behavior of the moments of the density of the blob along a Lagrangian trajectory 
is determined by the rate function H{a). Of course, 70 = 0. Eq. (12. 8p implies, in turn, that 7_i = 
since the minimum of H'{~a) is equal to zero. In the incompressible flow, all jq vanish. 

There are other quantities that may be expressed via the rate function H{a) of the large deviations 
of the stretching exponents including: 



-"^At long distances and not too long times, the effect of molecular diffusivity my be disregarded in practical situations. 
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rates of decay of moments of the tracer in the presence of diffusivity 
multifractal dimensions of the natural attractor measure [501 116] 

evolution of the moments of the polymer stretching B and the onset of drag reduction in polymer 
solutions El [261 



Conclusion. In smooth dynamical systems cumulation of local stretchings, contractions and rotations 
leads to the linear growth of the stretching exponents with the asymptotic rates equal to the Lyapunov 
exponents. If the Lyapunov exponents are all different then one expects that the statistics of the stretching 
exponents exhibits at long but finite times a large deviations regime characterized by the rate function 
11(3). Such a function contains more information about the short-distance long-time properties of the flow 
than the Lyapunov exponents and it enters the determination of several asymptotic transport properties 
of turbulent velocities at intermediate Reynolds numbers. The rate functions H for the direct and time- 
reversed flows are related by the Gallavotti-Cohen type symmetry \2.8\) . 
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3 Particles in smooth Kraichnan velocities 



In 1968 Robert H. Kraichnan initiated the study of turbulent transport in a Gaussian ensemble of velocities 
decorrelated in time. The statistics of such velocities is determined by the mean {v^(t,r)), that we shall 
take equal to zero, and by the covariance of the form 

{v\t,r)v'(t\r')) = 5[t~t')D''{r,r'). 

To imitate turbulent flows far from boundaries, one may assume: 

• homogeneity D^'-'{r,r') = D^-'{r — r'), 

• isotropy D^^r) = S'^ Di{\r\) + rV^ Dadr]), 

[ O(r^) for |r| <IC 77, 



scaling D*^ (r) - D'^ (0) 



0(|r|«) for r/ < |r| < L. 



The range where D^^ (r) — D*^ (0) is approximately quadratic is called the "Batchelor regime" and is 
used to mimic flows at intermediate Reynolds numbers, dominated by viscous effects, with the spatial 
velocity increments approximately linear in the point separation. The range where D^-' (r) ~ Z)'-'(O) 
scales with power ^ between and 2 mimics the developed "inertial range" turbulence at high Reynolds 
numbers, where the spatial velocity increments scale like a fractional power of the distance between the 
points (the Kolmogorov scaling is rendered by ^ = 4/3). Although not very realistic, the Kraichnan 
model of turbulent velocities has provided a very useful playground for theoretical and numerical studies 
of transport properties in random flows. It allowed to test old ideas about turbulent transport and, even 
more importantly, to come up with important new ideas [39) . 

3.1 Lagrangian particles and tangent flow in Kraichnan model 

In the Kraichnan ensemble, velocities are white noise (so distributional) in time and the Lagrangian 
trajectories are defined by the stochastic differential equation (SDE) 

dR = v{t,R) dt 

that we shall interpret with the Ito convention [691 164) . i.e. as equivalent to the integral equation 



R{t) = R{to) + lim f v{s,R{U))ds, 

where the limit of the Riemann sums defines the Ito stochastic integral Jl^^v{s, R(s)) ds. Note that the 
average value of the Ito integral vanishes since {v{s, R{ti)) = for s > ti. The Stratonovich convention 
would use the Riemann sums with v{s, R{ti)) replaced by v{s,R{-^(ti + ti+i))) but it gives the same 
result if diD^-' {r,r) = 0, which holds e.g. for homogeneous and isotropic velocities. 

For a regular function /, the Ito stochastic calculus gives the SDE for the composed process with an 
additional (Ito) term as compared to the standard chain rule: 

dfiR) = v{t,R)-Vf{R)dt + ^D"{0)\/,\/.if{R)dt 

with the assumption of homogeneity. For the averages, one obtains the differential equation: 



dt 

which may be rewritten in the form 



^(f{R)) = {^D^'m.V,f{R)) 



d_ 

di 



f{r){S{r-R{t;to,ro)))dr = J ^D'' {0) f{r) {S{r - R(t;to,ro))) 



dr 



or, stripping the last identity from arbitrary functions /, as 

^^{5{r-R{t;to,ro))) = ^D'' {0)\7 ,,\7 {S{r ~ R{t;to,ro))) . 

With the additional assumption of isotropy, iD*^ (0) = DqS''-' , and we infer that the probability to find 
the Lagrangian trajectory at time t at point r is that of the diffusing particle: 

P(to,ro;t,r) = {5{r ^R{t; to, ro))) = el*-*°l^°^ Vo, r) = (4^jj^|/_ ^p|)d/2 ^""^^^ ■ 
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The constant Do is called the "eddy diffusivity" . One infers that in mean, the Kraichnan turbulence 
causes difltusion. In particular, for the mean density and the mean scalar. 



{n{t,r)) = {Jn{to,ro)S{r~R{t;to,ro))dro)= y n(to, ro) e^'-'^'-^"^ (m, r) dro , 
{0{t,r)) = { f d{to,ro) 5{ro-R{to;t,r)) dro) = / ^(to, ro) e<*-*«'°"^ V> n)) dro 



for t > to, see Eqs. and (|1.10p . An initial blob of scalar or density diffuses in mean as watched in 

the laboratory frame. Of coarse, the evolution of the average density does not contain information about 
the dynamical behavior of the density fluctuations. A way to catch a glimpse of the latter is to look at 
the evolution of the blob from the point of view that moves with the fluid, as was discussed in Sect . 12.31 

To this end, let us study the dynamics of the small (infinitesimal) separation between the Lagrangian 
trajectories that satisfies the SDE 

dSR = {5R-'Vv){t,R{t))dt = j:{t)5Rdt, (3.1) 

where S%(t) = {'VjV^){t,R{t)). Again, we shall consider this equation with the Ito convention but the 
Stratonovich convention would give again the same result under the assumption that diD^-' {r ,r) = 0. 
For a function of &R , the Ito calculus gives: 

df{5R) = {5R-Vv){t,R{t))-^f{5R)dt- i<5i?*(57?' (VfcViD'^)(0) (V.Vj/)(<5Ji) di 

in the homogeneous case and for its average, 

^(/(5fl)) = {\5R''5R' C^, {VN,i){5R)) , (3.2) 

where C^^ = — VfeVjZ)*-' (0). The last result is the same as if SR{t) satisfied the linear Ito SDE 

dSR = S(t)5Rdt (3.3) 

with the matrix-valued white noise S{t) such that {S{t)} = and {S^k{t) S-' i(t')) = C^^ &{t — t'). Here 
taking the Ito prescription is essential since although the Ito and Stratonovich prescriptions give the same 
result for Eg. (|3.1[l . they do not agree for Eg. (|3.3p with S''j{t) being the white noise! Similarly, the 
statistics of the tangent process W{t;ro) for a fixed initial point vq may be obtained by solving the 
linear Ito SDE 

dW = S{t)Wdt (3.4) 

with the white noise S{t) and W{0) — Id. The decoupling of the Lagrangian trajectories R{t;ro)) from 
these equations is the great simplification of the Kraichnan model! 

The further simplification is the decoupling of the natural measures from the statistics of W{t; v\v)) 
for the homogeneous Kraichnan flows on a periodic box V. Due to the homogeneity, such measures 
n{dr\v) have to satisfy the relation 

{n{dr\v)) = 

Now, since n{dr\v) depends on the past velocities and, for t > 0, W{t,r\v) depends on the future ones 
(independent of the past ones in the Kraichnan ensemble), we infer that 

{ j^f{W{t-r\v)) n{dr\v)) = ^ {f{W{t;r\v))) {n{dr\v)) = {f{W{t-r\v)) ^) = {f{W{t-ro\v))) . 

It follows that the statistics of W(t\r\v) for t>Q (and hence also of p{t;r\v)) with {r,v) distributed 
with the natural invariant measure A'^(dr,d?;) coincides with that of the solution of Eq. (|3.4p with the 
white noise S{t) and VK(0) = Id (and of the corresponding p(t) ). In other words, 

Pt{S) = Pt{a) for t > 0. 

Since the Kraichnan ensemble is time reversible, we also infer from Eq. (|2.10|l that 

Pt{~5) = Pt{a) e'2:ti-. 

for all t > 0. 
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3.2 Multiplicative large deviations in Kraichnan velocities 
3.2.1 Isotropic case 

The specific form of Pt{S) depends on the covariance C^-J reflecting the single-point statistics of the 
velocity gradients. In the isotropic case, 

with the right hand side specified, up to normalization, by the "compressibility degree" 

In general, < p < 1. Vanishing p corresponds to incompressible velocities, whereas p = 1 to gradient 
ones. The normalization of C]^^ is set by the time scale that we shall define by the relation 

2dm _ 2d _ 2 ^gg^ 



{(V,t;J)2) Cf^ 213 + -yd' 

It may be thought of as mimicking the Kolmogorov time scale in the Kraichnan fiow. By the Ito calculus, 
for W{t) solving the SDE (|3.4p . one obtains the relation 



iim)) = {wrw^cz^^jm) . (3.7) 



Any function f[p) of the stretching exponents may be viewed as a function / on the matrix group 
GL{d) such that 

fiO'WO) = f{W) for 0',0(£0{d). 

To calculate ^(/(p)), it is enough to compute the operator inside the average on the right hand side of 
Eq. (|3.7p in the action on such functions. A straightforward although tedious calculation gives [19) : 

^ (fip)) = {^fip}) 

for the generator 

After the conjugation by the function .F(p) = e~ 2 Pi {Yli<_j sinh(pi — Pj))^^^, the operator Jiff becomes 
the Calogero-Sutherland-Moser Hamiltonian of d quantum particles on the line [65) : 



const. = -.J^fcsM 



Since, with the initial value p(0) — 0, 

(/(pW)) = / m) Pt{s) da, 

it follows that 

Pt{S) = te'-^(6,ta) = t lim JP(po)"' e-'-^^SM (^^^ j-) Tita) . 
Po-»o 

The spectral representation of the heat kernel e"*"^'^" is known explicitly [65) . Its saddle point calculation 
gives the large deviation form of Pt{p). The latter may be also found directly [7] since 

Pt(a) ~ te'^°°(0,ta) = const, e"'-^'"^' , 

where Jffoo is the second order differential operator with constant coefficients obtained from ^ by 
replacing cosh(pi — pj) by ±1 for i^j, respectively. An easy calculation gives then: 

with the Lyapunov exponents 

. _ (d+p(d-2))(d-2i + l) 



(d- l)(d + 2) 



P 
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equally spaced. In particular, Ai > for p < d/A (the chaotic phase) and Ai < for p > d/4. The 
change of sign of Ai induces a change in the transport properties (direct cascade of the the scalar for 
Ai > versus inverse cascade for Ai < The rate function H satisfies the Gallavotti-Cohen 

relation (|2.8p with H' = H, as may be easily checked. 

The exponents 7^ of Eq. (|2.17|l determining the asymptotic behavior of the density fluctuations in 
the Lagrangian frame are easily computed to be 

^qTri = pdq[q + l) . 

Note that 7^ > for q > Q and p > so that in the compressible homogeneous isotropic Kraichnan 
flow the positive moments of the density grow exponentially in the Lagrangian frame. 

3.2.2 Non-isotropic case 

The equal spacing of the Lyapunov exponents and the Gaussian form of the large deviations of the 
stretching rates is particular for the Kraichnan flow that is homogeneous and isotropic at short distances. 
On a periodic box, this is not a generic situation. For example, on a periodic square, typically, the 
covariance of the matrix-valued white noise S{t) will only have symmetries of the square and will be 
given by the expression 

= 2a5^ + + 5l5i) + ^S^^S.i , (3.8) 

where is equal to 1 if i — j — k — I and vanishes otherwise. The compressibility degree (|3.5|l and 
the Kolmogorov time p.6|) are now given by 

_ 2q + 3/j + 7 _ 1 

^ ~ 2(a + /3 + 7) ' " a + /3 + 7 

and one may introduce the "anisotropy degree" 

omega — |a|r,, . 

The Lyapunov exponents may be expressed in this case by the complete elliptic integrals [49] : 

3-2p + ^ E{k) 3^2p + uj E{k) 
^^"^ = " ^ + 2 Kik)' ^^r, = 1 - 2p , 

where k^ — , depicted below as functions of < < 1 for three values of the compressibility 

degree p = 0, p = ^ and p = 1 (the middle green lines give + A2)r,, — —p). 




If Ai > A2 then the large deviations of the stretching exponents are controlled by the rate functional) 
H{ai,(J2) = ir,,p~^(CTi + CT2 + 2pr~^)^ 

+ max[i/(ai-(72)-(l+a;)r,7V(// + l) + (3-2p + a;)r~^£;,,o(fc^)], (3.9) 
where Ei,fi{k'^) is the groundstate energy of the Lame-Hermite one-dimensional Schrodinger operator 

- 4-^ + u{u + l)k''sn^{u,k) (3.10) 
du^ 

acting on periodic functions of u with period 2K{k). sn(u, fc) is the Jacobian sine-amplitude function 
|49] . If Ai = A2, which happens for lu — 1 = p, then 

H{ai,a2) = ^-r,,((7i + (72 + 2r~^)^ + ir,,(cri - 0-2)^ (3.11) 

and cannot be obtained as the limit of the previous formula (the time-scales at which the large deviation 
regime sets in diverge as (Ai — A2)t,, —> [31]). Except for the last case, the large deviations of (pi — P2) 
are non-Gaussian in the presence of anisotropy, i.e. for uj > 0. 
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3.3 Inertial particles in Kraichnan model 

The inertial particles dynamics in the Kraichnan velocities is given for the vanishing diffusivity k by the 
SDE 

dR^Udt, dU ^ ^{-U + v{t,R))dt, (3.12) 

see Eq. (|1.4|l . Both Ito or Stratonovich conventions give for it the same result. Small phase-space separa- 
tions between two inertial particles satisfy, in turn, the equations 

dSR = SVdt, dSU = ^{-SU + T,{t) 5R) dt , (3.13) 

where, as before, S'j(t) = V jV^{t, R{t)). The matrix W{t;r,u) propagating {5R,5U) along the tra- 
jectory starting at (r,u) satisfies then the SDE 

dW = (i^^^^ ^'{)wdt. (3.14) 

As before, if the velocity statistics is homogeneous then E(f) may be replaced by the matrix- valued white 
noise S{t) as long as we are interested in the statistics of W{t; r, u) for fixed (r, u). The same is true if 
we average over the latter points with the natural measures n(dr,du) since the latter average decouples 
due to the independence of the past and future velocities). For the homogeneous and isotropic Kraichnan 
velocities, on a periodic box V, the ensemble average of the natural measures is equal to the Majcwell 
distribution: 

(n(dr,du)) = ^'^^'^.^ — e ^do^ dr dit . 

Note that Eq. 1)3. 14[) implies that det W(t) = eT'^^^^ so that the corresponding stretching rates o"i(t), i = 
1, . . . , 2d, satisfy the relation 

2d 

^a,(t) = -^. (3.15) 

1=1 

A simple argument [40] shows that the Gallavotti-Cohen type symmetry (|2.8p for the large deviations 
rate functions reduces now to the identity 

//(r-'l-a) = H{3) (3.16) 

with 1 = (1, . . . , 1) holding for a satisfying the relation (|3.15|l . 



3.3.1 One-dimensional case 

For the one-dimensional Kraichnan velocities, setting (5r(t) — ip(t) e*^^'^'^^ permits to rewrite Eqs. (|3.13|l 
in the second order form [80] 

-i, + V{t)f = (3.17) 

with V{t) = iE(f). Again, V{t) may be replaced by a white noise with the covariance 

{V{t)V{t') = ^S{t^t') (3.18) 

Upon viewing i as a space variable, Eq. (|3.17p becomes the one-dimensional stationary Schrodinger 
equation in the delta-correlated potential V{t) at the energy E — — j^r, a well known model for the 
one-dimensional Anderson localization. This is an old problem whose solution goes back to [51j . see also 
|59) . One considers the process x{t) — tp{t)/^f}{t) that satisfies the SDE 

dx = -(a;^ + E- V{t)) dt . (3.19) 

Its trajectories reach — oo in a finite time but reappear immediately at -|-oo . Such jumps correspond 
to a passage of ^{t) through zero with ^ or to a crossing of close inertial particles in space. The 
jumping process possesses the stationary probability measure 

X 

n{dx) = j e-^-^(h'+Ey)Ay)dx, (3.20) 
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with the density ~ x~'^ for large ( Z is the normahzation constant). The locahzation Lyapunov 
exponent A, equal to the mean value of x{t) = ^ ln(|?/'(f)|), is always positive signaling the permanent 
localization in one dimension. The Lyapunov exponents for the inertial particles are -^i = ^ ~ 2t ''■'^'^ 
A2 = —A— ^ (recall the relation between ^(t) and 5r(t)). A can be expressed [591 by the Airy functions 
[49] and one obtains: 

> -r A.M , , Ai'(^^)Ai(^^) +Bi'(^^)Bi(^.^) N 

Ai.^r, = 4n - n ± Ai(n^) + Bi(^^) ) ' ^^'^^^ 

for u = l;St~^^^ and the Stokes number St — —. For large St, Ai is positive and decreases as St~^^^. 
For St small, it is negative and approaches when St the value —1 for the Lagrangian particles. 



3.3.2 Two-dimensional case 

For the two-dimensional homogeneous and isotropic Kraichnan velocities, one may treat 5r (t) e*/^'^^' = 
ip(t) as a complex-valued process which still satisfies Eq. (|3.17|) with V{t) a complex valued white noise 
with the covariance [671 162] 

{V{t) V{t')) = 5{t - t') , {V{t) nF)) = ^5{t^t'). (3.22) 

For the complex-valued process z{t) = ip{t)/xp{t) one obtains now the SDE 

Az = -{z' +E- V{t)) dt (3.23) 

which may be shown to possess an invariant probability measure, this time without an explicit analytic 
form but decaying like for large \z\. The top Lyapunov exponent Ai is equal to the expectation 

of Kez in this measure. It was studied numerically in [621 1521 [H] . It seems again to behave as St~^^^ 
for large Stokes numbers. The reference [14] contains also the numerical results about large deviations of 
Rez{t). 



Conclusion. In smooth Kraichnan velocities with symmetry properties, the tangent process for La- 
grangian particles is related to integrable quantum-mechanical systems. This relation gives an analytic 
access to the large deviation rate functions of the stretching exponents. For inertial particles, similar 
relations to one- dimensional models of the Anderson localization permit to find analytically the Lyapunov 
exponents in one-dimensional flows and facilitate their analysis in two (and three) dimensions. 
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4 Particles and fields in rough Kraichnan velocities 

The previous discussion pertained to the transport of particles in the Batchelor regime with smooth ve- 
locities and small particle separations. We shall pass now to the study of particle dynamics at separations 
in the inertial range of scales where velocities become rough. We shall concentrate on the motion of 
Lagrangian particles. For recent ideas and numerics about the motion of inertial particles in rough flows, 
see [Smi]. 

Let us start by considering the simultaneous motion of A'^ Lagrangian particles Rn{t) in the homo- 
geneous Kraichnan velocities. A function f(R) of their joint positions R{t) = {Ri{t), . . . , i?jv(t)) evolves 
accordingly to the Ito stochastic equation 

N N 

df{R) = ^«'(t,flOVBj^/(^)df + ^ ^ D'^{R^-R„)Vj,,^V^,^f{R)dt. (4.1) 



2 

m.n — 1 



For the expectations, this gives: 

^^{fm) = iiMMfm)). (4.2) 

4.1 Particle dispersion 

First, let us study the evolution of a function of the separation AR{t) = Ri{t) — R2{t) between two 
trajectories: 

^(/(Aii)) = {{D^^iO) - D'^{AR))\7,\7,fiAR)) . (4.3) 

In particular, we shall look at functions of the inter-particle distance A = |Ai?|, called the "particle 
dispersion" . In the Batchelor regime, 

D'' (0) - D'' (AR) ^ ^AR''AR'VkViD''{0) 

and Eq. (|4.3p reduces to the same equation as for {^f{5R)), see the relation p.2|) . For functions of the 
particle dispersion A = | AJi|, one obtains after a short calculation, assuming also isotropy of the velocity 
ensemble: 

where Ai = ^0^7 — /3 is the first Lyapunov exponent. For the PDF of the distance A this gives 



i.e. a log-normal distribution with Ao = |Ai?(0)j. Note that 



{5{A - \AR{t)\)) = ^ ^ - e 2o+7)t ^ ^ 



^liin^ {5{A-\ARm)=5{A) 

so that when the initial distance Ao of the trajectories tends to zero, so does their time t distance. The 
latter behavior is imposed by the uniqueness of the Lagrangian trajectories in each velocity realization, 
given their initial point (and assuming that the molecular diffusion is absent) . It is consistent with both 
exponential separation (when Ai > 0) and exponential contraction (when Ai < 0) of close trajectories. 

Now, let us consider the Kraichnan model in the inertial range. Here D'^ (0) - D'^ (r) = 0(r«) and 
the isotropy imposes the form of the tensor: 

D''{0) - D''{r) = Pr'r'\r\^-'' + ^-jS'^r^ 

in the way fixed, up to normalization, by the compressibility degree 



= lim 



^0 VjVjL»"(r) 2/3 + d-/ 



between and 1 (taking ^ = 2 reproduces the smooth case formula (|3.5p ). The equation for the evolution 
of the two-particle dispersion A takes now the form 
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for = ^W+V' ^ T+k~'^- The distance process A{t) is the one-dimensional diffusion with the generator 
Ala and the transition PDF 

(5(A-A(t;Ao))) = e*"»(Ao,A). 

It is instructive to change variables from A to a; = A^"*''^ since the process x{t) is the Bessel diffusion 
with the generator proportional to x^~^''^'^dxX^''"~^dx, i.e. to the radial Laplacian in Deff dimensions, 
where 

2(a + l-g) 

In other words, a;(t) behaves as the norm of the Brownian motion in (continuous) dimension D^s- 
The short-distance behavior of the Brownian motion below and above two dimensions is different. In 
particular, for < DcH < 2, i.e. for 



1 d-2 1 d 
Pc = + o<P<75-=Pi 



2^ 2 ° e 



2 



the generator of the Bessel process admits different boundary conditions at a; = 0. They correspond to 
different behaviors of Lagrangian particles when they collide. 

4.2 Phases of Lagrangian flow 

There are three different phases with trajectory behavior and different transport properties in the inertial 
range, depending on the level of compressibility [481 1331 [57] . 

4.2.1 Weakly compressible phase 

For < p < Pci i.e. for D^s > 2, there is only one possible boundary condition at zero: 

{A''^/)(0) = 

and it renders Ma self-adjoint on the Hilbert space L^(R4-, A''~^dA). The process A{t) has a reflecting 
behavior at zero and 



lim (5(A - Aft; Ao))) = const. 

An^O 



t 2-5 



The formula says that the distance between two particles at time t is smoothly distributed even if their 
initial distance is zero: the trajectories separate in finite time (unlike in the chaotic regime with exponen- 
tial separation where close particles take longer and longer to attain a sizable separation)! This implies 
the "spontaneous randomness" of the Lagrangian flow TS^: the sharp distribution 5{r — R{t;to,ro\v) 
of the time t position of the trajectory starting at a fixed point in a fixed velocity field realization is 
replaced by a diffused distribution P{to,ro;t,r\v) that is not concentrated at one point, as illustrated 
below: 




This is possible because for D{0) — D(r) cc |r|* the typical velocity realizations v{t,r) are only Holder 
continuous in space, with exponent smaller than ^/2. In such velocity fields, the solutions of the trajectory 
equation ^ — v{t,R) in a fixed v are not uniquely determined by the initial condition but form 
a "generalized fiow": a stochastic (Markov) process with transition probabilities P{to,ro;t,r\v). The 
latter have been constructed rigorously by Le Jan and Raimond in [57) . 
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The spontaneous randomness of the Lagrangian flow implies the persistence of dissipation in the scalar 
transport. Note that for the incompressible flow where 



P(t,r;to,ro\v) dro = J P(t,r;to,ro\v) dr = 1 
and the unforced scalar evolves according to the equation 

0{t,r) ^ J P{t,r;to,To) e{to,To) dro, (4.4) 

the inequality 

0< J P{t,r;to,ro\v) {9{to, ro) - 9{t,r)f dro dr = J e(to,rof dro - j e{t,rf dr 

implies that the norm of 9 (the scalar "energy") cannot grow. Besides, it is conserved if and only if, 
for each r, 6{to, ro) ~ 0{t, r) on the support of P{t, r; to, -Id), so if and only if the flow is deterministic. 
The persistent dissipation of the scalar energy in generalized flows is a short-distance phenomenon. It 
leads to the direct scalar-energy cascade where the scalar energy supplied at long distances is transfered 
without loss towards shorter and shorter scales and finally dissipated on the infinitesimal ones (in the 
limit of vanishing molecular diffusivity k) |55I 139) . 



4.2.2 Strongly compressible phase 

For p > Pc, i-e. for DeH < 0, the only possible boundary condition for Ma is /(O) = and it again 
renders Ma self-adjoint in L'^(K+, A'^~^dA). Now, A{t) is absorbed at zero and 

((5(A - A(t; Ao))) = regular 4- const. (5(A) , 

see [48]. In this phase there is a positive probability for the Lagrangian trajectories that start at positive 
distance Ao to collapse together by time t. When Ao — > then the contact term dominates: 

^lirn^ <5(A-A(t,Ao))) = S{A) 

which signals that the trajectories are deterministic: 

P{to,ro;t,r\v) = 5{r - R{t;to,ro\v)) 
for single trajectories R{t;to,ro\v) collapsing for different ro at later times, as on the illustration: 



R=R' 





This is again a non-standard behavior. Since trajectories are unique, there is no persistent dissipation of 
the scalar energy and, in the presence of a steady source, the latter exhibits an inverse cascade (towards 
long distances) [48] . 

4.2.3 Intermediate compressibility phase 

For pi < p < Pc, i.e. < -Dcff < 2, both refiecting and absorbing boundary conditions are possible, 
the first one is selected by adding very small diffusivity, the second one by considering velocities smeared 
at very small distances (small viscosity) [33]. The other permitted boundary conditions are the "sticky" 
ones [i5] : 



A'^^/KO) 
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for /i the amount of "glue", < fj. < oo. These conditions render Ma self-adjoint on the space 
L^(R+, (A""^ + fi5{A))dA). The particles stick together and separate immediately, never spending a 
finite interval of time together. Nevertheless, they spend a finite portion of the total time glued together. 
The persistent dissipation is still present and leads to a direct scalar-energy cascade of the forced scalar. It 
induces however an anomalous scaling of the 2-point scalar structure function {{9{t, r) — 6{t, 0))^) cx r^~" 
(instead of the dimensional result oc r^~*, see below). The sticky condition corresponds to fine-tuned 
small diffusivity and small viscosity [SH |46] . 

For pi < p < Pc, the transition probabilities P{to,ro;t,r\v) were constructed rigorously by Le Jan 
and Raimond |57l 158) . for the reflecting boundary condition, the absorbing one and, in one dimension, also 
for the sticky ones. In the first case, they depend on v only, whereas in the second and third cases, on the 
(white noise) v and on the additional "black" (non-standard) noise [77 j that decides what the particles do 
when they meet. Constructing the transition probabilities of the generalized flow for the sticky boundary 
conditions in more than one dimension is still an open problem. 



4.3 Zero-mode mechanism for scalar intermittency 

The statistics of fully turbulent velocities exhibits the phenomenon of intermittency, characterized by 
frequent appearance of strong signals interwoven with weak ones both in time and in space domain [42j . 
Similar, but even more pronounced effects characterize the statistics of advected scalars [31 [B^- One of 
the motivations behind studying scalar advection in Kraichnan velocities, whose Gaussian statistics does 
not exhibit any intermittency, was a search for a mechanism of this enhanced scalar intermittency [56[ . 
We shall briefly describe such a mechanism discovered subsequently [281 1471 174) and its interpretation in 
terms of the related flow of Lagrangian particles. 

4.3.1 Passive scalar in Kraichnan velocities 

For the scalar advected by smooth velocities, 

9{t,R{t;to,ro)) = 9{to,ro) = const. 
In the Kraichnan velocity ensemble this means that the scalar field has to satisfy the Ito stochastic PDE 

dt9{t,r) = -{v{t,r)dt) ■'V9{t,r) + ^DoV^9{t,r) dt . 
Indeed, the Ito rule gives then 



dt9{t,R{t)) ^- {v{t,R{t))dt) ■V9{t,R{t)) + ^DoV^9{t, R{t)) dt 
+ {v{t,R{t))dt) ■'V9{t,R(t)) + ^Do'V^9(t,R{t))dt 
- DoV^9{t,R{t))dt = 

(the last term in the intermediate equality came from the mixed increments of 9{t) and of R{t)). The 
corresponding Stratonovich stochastic PDE would be 

dt9{t,r) = ~{v{t,r)dt) ■V9{t,r) 

(without the diffusion term). Similarly, for the scalar with a source, one obtains the Ito equation 

dt9{t,r) = -{v{t,r)dt) ■V9{t,r) + ^Do'V^9{t,r) dt + g{t,r)dt (4.5) 

Let us suppose that g{t,r) is also a Gaussian process, independent from velocities, with mean zero and 
covariance 

{g{t,r)g{t',r')) = 5(1 ~ t') xir - r') . 
For the product n!^=i ^(^i ^'Oi '^^^ obtains then from Eq. ([4.5|l the Ito SDE 

JV JV JV 

dtl[9(t,r^) = (i;(f,r„)df) ■ Ve(t,r„) JJ 9(t,r^,) + ^DoV^9{t,r„) JJ 9{t,r„,)dt 

n—l n — l n'j^n n — l n' 

-f ^ D''(r^ -r„)V«6l(t,r,„)Vj6l(f,r,0 J]^ 9{t,r,,,)dt 

N 

+ ^((7(f,r„)df) ]^ 6l(f,r„0 + ^ x(T-n-»-n') H ^(i=^n')dt. 

n — l n' n<.n' n" ^n^n' 
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For the expectation values, this gives the relation 

JV JV 
n — l n — 1 n<n' n" ^n^n' 

where the operator is the same as in Eq. (|4.2p for the evolution of expectations of functions of 

positions of A*' Lagrangian particles. 



4.3.2 Zero modes and anomalous scaling 

Under steady forcing, the scalar statistics reaches a stationary state in the weakly compressible phase 
with direct scalar energy cascade. In this state. 



The above equations determine the scalar correlations ( Y\ ^(^i '''■",)) inductively in A'', up to zero modes 
of the operators . It appears that if the source g is concentrated mostly at long distances then the 
short distance scaling of the scalar "structure functions" S^{r) = {{0{t,r) — 9{t,0))^'^, 



(4.6) 



S„{r) oc |r|'^" for small jr| , 

is determined by the contribution of the scaling zero-modes of A4^ : 

M^f^^O, ./■„(A^) = A««/„(fl). 

Recall that the functions of the TV Lagrangian trajectories evolve by Eq. (|4.2p . The zero-modes of A4^ 
correspond to functions f{R{t)) of the TV-particle process that are conserved in mean: 

{fmt^Lo))) = ./(rio). 

They are martingales of the effective multi-particle diffusions. 

The scaling zero-modes of A4j^ giving the leading contributions to the structure function Sj^{r) have 
been found for homogeneous and isotropic flows in the leading order of the perturbation expansion in 
powers of 5 [171 [18] and of i [28l[27] in the years 1995-6. The result was: 



K(9 - n - iV(iV-2)(l+2p); 
2 'i' 2(d+2) 

iV(iV-2)(l+2p) i 
2d 



f(2-C)- 



-o{e) 



(4.7) 



anomaly 



Terms to the 3^^ order are known [T]. The fact that < -^(^2 means that 



&(r)^/2 



becomes large for small jrj. This signals the presence of long tails in the stationary distribution of the 
scalar increments (9{t,r) — 9{t,0)), getting more and more pronounced at short distances. This is a 
signal of small scale intermittency of the scalar statistics observed also in other turbulent states. The 
perturbative results (|4.7p were confirmed by numerical simulations, see the figure from [43] 
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representing the results for 2C,2 — ^4 as a function of ^ in two- (upper points) and three-dimensional (lower 
points) incompressible Kraichnan model, with the broken line in the lower left corner representing the 
perturbative result (|4.7p for d = 3. Zero-modes were also shown numerically to be also responsible for 
the intermittency of scalar transport by the inverse cascade of 2D turbulence [25) . 

Conclusion. In the non-smooth Kraichnan velocities modeling the inertial range turbulence, Lagrangian 
particles exhibit non-standard behaviors impossible in differentiable dynamical systems but expected to 
be common in non- differentiate ones with Holder- continuous vector fields. These behaviors include ex- 
plosive separation of trajectories generating their spontaneous randomness and an implosive collapse of 
deterministic trajectories, with possible stickiness. Such unusual behaviors of trajectories are source of 
non- equilibrium transport phenomena involving cascades with non-zero fluxes of conserved quantities. 
They condition the appearance of short-scale intermittency tn the direct cascade of scalars advected by 
fully turbulent flows. The Kraichnan model allowed to associate the scalar intermittency to hidden statis- 
tical conservation laws of multi-particle evolution: the "zero modes". The zero-mode scenario permitted 
perturbative calculations of the anomalous scaling exponents of the scalar structure functions in this model. 
Similar mechanism is believed to be responsible for intermittency is other problems with dynamics governed 
by linear evolution equations y^. 
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5 Final remarks 



The passive transport in turbulent flows is related to the behavior of Lagrangian or inertial particles 
carried by the fluid. In the regime where velocities are smooth in space, i.e. for intermediate Reynolds 
numbers, the particle dynamics provides examples of random differential dynamical systems and may be 
studied with tools of the chaotic dynamical systems theory. It requires, nevertheless, subtle information 
about rare dynamical events (multiplicative large deviations). The statistics of rare fluctuations obeys 
fluctuation relations, sharing this property with other non-equilibrium systems [3U] . 

In the regime of large Reynolds numbers, i.e. in the inertial range, the motion of particles may 
be viewed as providing an example of random non-differentiable dynamical system. As the work on 
the Kraichnan model has shown, such systems lead to unconventional flows with explosive trajectory 
separation or implosive trajectory trapping. The exotic behavior of particle trajectories has dramatic 
effect on the transport properties and conditions the appearance of cascades of conserved quantities and 
of intermittency related to hidden conservation laws. 

Some of the lessons from studying passive turbulent transport remain valid for reactive particles like 
water droplets or chemical or biological agents |38l 176) . How many of them may be be transformed to 
insights about turbulence itself remains to be seen [2]. 
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